ADP-94-27/T166 



December 1994 



A Study of Structure Functions for the Bag Beyond Leading Order 

P.M. Steffens and A.W. Thomas 



Department of Physics and Mathematical Physics 
University of Adelaide 
Adelaide, S.A. 5005, Australia 

Abstract 

There has recently been surprising progress in understanding the 
spin and flavor dependence of deep inelastic structure functions in 
terms of the same physics needed in the simple quark models used for 
hadronic spectroscopy. However, the corresponding scale is usually 
very low, casting doubt on the use of leading order QCD evolution. 
We show that the conclusions are not significantly altered if one goes 
to next-to-leading order. In particular, the excellent agreement with 
unpolarized and polarized valence quark distributions is retained. 
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1 Introduction 



Deep inelastic scattering continues to provide a wealth of surprising and 
challenging information concerning the structure of nucleons and nuclei. The 
nuclear EMC effect, now more than a decade old provided a dramatic 
challenge to our conventional view of nuclear structure 0. In particular, it 
suggested a systematic change in the valence quark structure of a nucleon 
inside matter. This clearly goes to the heart of our understanting of nuclear 
structure. More recently the European Muon Collaboration revealed |0] an 
unexpected deviation from the Ellis- Jaffe sum rule [Q, which became known 
as the "spin crisis" . Finally the New Muon Collaboration has confirmed 
suggestions of a violation of SU(2) flavour symmetry in the nucleon sea. 

All of these observations demand theoretical interpretation. In the ab- 
sence of reliable calculations using lattice QCD, and with QCD sum rules 
yielding only a few moments, one is naturally led to the problem of relat- 
ing low energy quark models to the parton distributions measured in deep- 
inelastic scattering. Given that deep inelastic data is valence dominated at 
low-Q^, it was suggested fairly early ^ that the natural connection was to use 
the quark model to calculate the leading twist parton distributions at some 
low scale, /i^, and to evolve them to a higher scale, Q^, using the renormal- 
ization group equations (RGE) and then compare the resulting distributions 
with data 0. In the context of the bag model, Jaffe [§] and Le Yaouanc et 
al. were the first to implement this idea, followed by Hughes fl^ and Bell 
ll|] and others [|l^ . 



A major problem in these calculations was the lack of correct support 
for the parton distributions. ( A quark distribution has to vanish outside the 
region < x < 1.) This disease was cured some time ago in a series of papers 
by Signal, Thomas and Schreiber |14| . The proposal was to construct the 
parton distributions using the formal analysis of Jaffe []T3[, and to ensure 



energy-momentum conservation before any approximations were made. (To 
better appreciate the importance of correct support in the distributions, the 



reader is referred to ref. |16|.) In ref. |14| a detailed study of the MIT quark 



distributions was made. Various distributions, polarized and unpolarized, 
were calculated at the model scale and then evolved in leading order (LO), 
using the RGE, to the scale of the experiment. A very good, qualitative 
agreement was obtained, which was rewarding and stimulating given the 
simplicity of the model. It became clear that high energy experiments could 
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14 



are not 
ref. 



be, to a large extent, explained using as input low energy physics. 

Of course, bag model calculations using the method of ref. 
the only ones that can be found in the literature (see, for example 

Moreover, there are calculations based on other models. For example, 
there are works on the colour-dieletric model [|18| or on the non-relativistic 
quark model ITB, O . However, besides the fact that some of the calculations 



suffer from the problem of poor support, it is known that the parton 
distributions constructed from Gaussian wave functions do not have the right 
behavior as x approaches 1. Recently, there have also been attempts to 
formulate the problem in terms of relativistic quark-nucleon vertex functions 
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parametrizations |2^, ^ 
some QCD sum rules 12 



Thus far, all analyses of the x dependence of parton distributions within 
quark models have been performed only in leading order. There are nu- 
merous studies of higher order evolution of unpolarized data in the form of 

or higher order corrections to the first moment of 
but no study of models involving both polarized 
and unpolarized data Besides the fact that corrections of higher order 
in the strong coupling constant ^^(Q^) are already important per se, we 
emphasise that a next-to-leading order (NLO) analysis is necessary because 
the scale of the quark model is generally very low so that the applicability 
of the leading order RGE is questionable. In this paper, we will apply the 
NLO evolution to the polarized and unpolarized valence distributions in the 
proton and to the x dependence of the Bjorken sum rule. We will compare 
the resulting distribution with the MRS |23| parametrizations and with pre- 
liminary data for the polarization of valence quarks in the proton [^. The 
paper is organized as follows. In Sec. II we briefly review the process of 
QCD evolution in NLO. In Sec. Ill we review the procedure for calculating 
quark distributions in the bag model. In Sec. IV we present some numerical 
results, while Sec. V is used to present our conclusions. 



2 Distributions at Next to Leading Order 

It is well known that one can write the moments of the structure functions 
as: 

^In ref.psl the x dependence of gip{x) was calculated in NLO. However, parametriza- 
tions of the data were used for the input 
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M„(g^) = EQ(W,^?)^n(/^')- (1) 

i 

The sum runs over all spin-n, twist-2 operators. In this work we will restrict 
ourselves to the nonsinglet sector because the complete, singlet anomalous 
dimensions in two loops for polarized scattering are not known. (Up to now, 
only the quark-quark and the quark-gluon part of the singlet Altarelli-Parisi 
splitting functions have been calculated p5| .) 



Equation (P is a direct result of the operator product expansion applied 
to the forward scattering of the photon from the hadronic state through the 
e.m. current, J^. The amplitudes A\{ii'^) are target dependent and involve 
non-perturbative QCD. They are related to the moments of the quark distri- 
bution in the target at the renormalization scale /i^. The Wilson coefficients, 
C^(Q^//i^, (yf), are target independent and calculable in perturbation theory. 
They carry all the information about the scale dependence of the structure 
functions and are well determined. The evolution in for the Wilson coef- 
ficients is given by the solution of the corresponding renormalization group 
equation which has the form: 



Here, giQ"^) is the running coupling constant (such that gi^Q"^ = /i^) = g), 
^^{g) is the anomalous dimension of the corresponding nonsinglet operator 
and P{g) is the QCD beta function. Expanding to second order they are 
expressed as: 

/'W = -''oi|^-A^. (3) 
,.,e'^).,<o.._|L (4) 

while the running coupling constant is determined by solving the transcen- 
dental equation: 



A2 Pot Pi 



I67H A 

Pot Pi 



(6) 
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The QCD scale parameter, A, is determined by comparing the theoretical 
calculations with the experimental data, and Po, Pi, ^C^)"-, '-)/(^)"- and are 
parameters determined in perturbation theory. Their form is listed in the 
appendix. Throughout this paper we will use A = 0.2 GeV (which is within 
the experimental errors For unpolarized scattering both the leading 

order (LO) p8| and next to leading order (NLO) p9[ coefficients are well 
known. One need only be careful with the renormalization scheme depen- 
dence: although Po, Pi and 7*^°^ do not depend on the scheme, 7^^)" and F" 
do, and both have to be calculated in the same scheme in order to obtain 



physically meaningful results BO]. For the nonsinglet operators, all coeffi- 



cients but are the same for both polarized and unpolarized scattering. 
The Wilson coefficient F„ was ffist calculated in the MS scheme by J. Ko- 



daira and collaborators 311. With this in mind we write the NLO evolution 



equation for the moments with the help of expressions (^) 



y(0)n 



/2/3o 



1 + 



167r2 
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(l)n 



(O)n' 



167r2 



Pi2_ 

2Po 2Pl 



•(7) 



This equation can be rewritten as: 



M„(Q2)=A„(g2)C„(l,^(Q2)), (8) 

which makes it clear that the moments of the quark distributions in NLO 
have the same dependence as in LO, but that this is not the same for 
the moments of the structure functions. In other words, y4„(Q^) gives the 
moments of the quark distributions and M„((5^) gives the moments of the 
structure functions. Expression also expresses clearly the scheme depen- 
dence of the quark distributions beyond leading order, which means that 
they are unphysical. From eq. (|^) is straightforward to verify that in NLO, 
QCD sum rules, like the Bjorken sum rule, pick up corrections. 



3 The Quark Distributions for the Bag Model 

The matrix elements 74„(/i^) can be written in terms of parton distributions, 
as showed by Jaffe |^2[|. This is done by defining the distribution as the 



5 



integral, in the light cone gauge, of the forward virtual quark-target scattering 
amplitude over all the parton momenta (in the light cone) but keeping the 
plus component fixed and equal to xp~^ (with x the Bjorken variable and 
p'^ the plus component of the target momentum). The parton distribution 
written in this way has support only for — 1 < x < 1. Explicitly, 

g(x)=p+^5(p+(l-x)-p+) I {n\^+mP) T, (9) 

n 

and 

r+oo 

An= / g(x)x"~Mx, (10) 



where Pn is the momentum carried by the spectators to the interaction (inter- 
mediate states). An equation very similar to equation (^ can be written for 
the antiquark distribution, but with "^"^ replacing \E'. There are two kind of 
contributions to (|TD|): one coming from the annihilation of one of the quarks 
in the nucleon and the other from the creation of a quark or antiquark in 
the nucleon. The first process is what we call the two quark contribution 
and it is the dominant one. We show this contribution within the model, in 
question, in detail below (see equ. (0)). The second contribution is a four 
quark one. As will become clear below, this kind of contribution will be used 
to fix the normalization problem of the quark distributions in the bag but it 
will not be calculated explicitly within the model. 

The approach followed in this paper to calculate q{x) will be that de- 
veloped by the group at the University of Adelaide [l^, [l^, in which the 
MIT ||33| bag model wave functions are used and the correct support for 
the distribution is fully assured. Basically, what is done is to project the 
moment vectors \p) in the coordinate space and calculate these overlaps us- 
ing the Peierls-Yoccoz aproach |3^. This ensures that the wave functions 



are momentum eigenstates, so attacking one of the major problems of the 
bag model, that is the lack of translational invariance. The Peierls-Yoccoz 
projection is regulated by the weight functions 0(p) given by: 



(p) |2= J dxe-'^-^ J dyij^{y-x)^{y) \ (11) 



Inserting a complete set of position states, converting the sum into an integral 
over the momentum of the intermediate states and using the projection given 
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above, we get for the two quark piece (full details can be found in ref. ||14|| ): 



^W4'(P.)P. (12) 



Tl M ^ 

X I n I d I n I 

y[Af2(l_:r)2-M2]/2A/(l-a;) 103(0) p 

Here is the spin- flavor part of the wave function of the initial state (at 
rest), Pf^rn makes the projection onto flavor / and spin projection m, M„ 
is the mass of the intermediate state and the Fourier transform of 
The expression for the four quark contribution has the same form, but 02 is 
replaced by 04. 

Once we have a model to calculate the matrix elements we can proceed 
to evalute the moments of the structure functions. Before doing so, a few 
remarks need to be made. The first deals with the spin of the intermediate 
states. In the case of a two-quark spectator system, the sum of the spins 
may be zero (a scalar diquark with a lower mass M^), or one (a vector 
diquark with a higher mass M„). The mass difference between these states 
is found to be around 200MeV in most quark models, and this is crucial 
to understanding the soin and flavour dependence of the distributions . 
Taking this into account, one can calculate the contributions to the quark 
distribution (|^) due to the scalar and vector diquark systems, by using the 
SU{6) wave function for The second issue is the form of the four quark 
contibution. It is known [|14| that the bag model has problems to fix the 
normalization of the quark distributions. Probably this is due the fact that 
we do not calculate the masses of the intermediate states in a way consistent 
with the model. To avoid this problem, we choose the four quark term 
to have the form (1 — xY and normalize it in a way that the overall quark 
distribution satisfies the normalization condition. Although this procedure is 
clearly phenomenological and not directly related to the model, we note that 
a term of this form will affect only the region x < 0.3. Furthermore its shape 
resembles very much the actual form of the the four quark term calculated 
using the same formalism as that used for the two quark contribution. 
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4 Results 



Now that we know how to calculate the matrix elements An we can proceed to 
apply the evolution equations to this particular model. The procedure we will 
follow is: first, we find a best fit for the unpolarized valence distribution by 
comparing the LO and NLO evolution against the MRS ||23| parametrization 
of the experimental data at 10 GeV^. In doing that we fix the set of model 
parameters, e.g., radius of the bag and masses of the intermediate states. 
Using the virial theorem in the bag, where the bag itself carries the same 
energy as each of the three quarks confined in the bag, we roughly expect 
M, ~ 750 MeV - 150 MeV and M„ ~ 750 MeV + 50 MeV. But this is, 
of course, only an approximation and we have some freedom around these 
values. We also determine the scale /z at which the model is supposed to be 
valid. Using this set of parameters, we can then calculate the predictions of 
the model for other quantities. Here, as we have only the nonsinglet evolution 
for NLO, we will restrict ourselves to the polarized valence distribution of 
the proton and to the difference between the spin structure functions of the 
proton and neutron. 

In LO it is found that for R = 0.8 /m, M, = 0.55 GeV, = 0.75 GeV 
and /i^ = 0.0676 GeV^, the unpolarized valence distribution of the proton 
in the bag model fits the MRS parametrization of the data. The result is 
displayed in Fig. |l] together with the bag model distribution at scale /i^. The 
agreement is excellent with no significant discrepancy. However, at such a 
low scale the coupling constant is q;((5^ = /i^) = 2.66. This is a rather large 
value (as found in ref. [0]), and raises doubts about the whole procedure. 
In NLO the parameters used to fit the unpolarized valence distribution have 
the following values: R = 0.8 /m, M, = 0.7 GeV, M„ = 0.9 GeV and 
/i^ = 0.115 GeV^. We notice that the radius for LO and NLO is the same 
but not the masses, being slightly bigger in the latter case. This is because 
in NLO one need not evolve as far, the scale /i^ is larger and so the masses 
of the intermediate states need to be larger (the parton distributions peak at 
a larger x for smaller intermediate state masses). In Fig. the NLO results 
are shown and again we see a very good agreement between the calculated 
valence distribution and the MRS parametrization of the data. What is 
remarkable now is the fact that the strong coupling constant drops to 0.77 
at /i^ = 0.115 in NLO. Although this value is still large, it is much smaller 
than the LO case, making the evolution more reliable. Both LO and NLO 
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moments are evolved to = 10 GeV^. 

Once we have determined the parameters of the model, we are in the 
position to calculate other quantities. In Fig. ^ we plot the polarized u va- 
lence distribution for the bag model in LO and NLO against the preliminary 



experimental data |26[. There is no real difference between the LO and NLO 
curves (besides, of course, the mass parameters). Due to the large errors in 
the data it is impossible to conclude anything for the region x < 0.1. But for 
X > 0.1 the theoretical curves certainly have the correct behaviour. More- 
over, the polarized valence distributions are pure nonsinglet and we might 
expect that the bag model would give, in this description of data as 

good as in the unpolarized case. This is because the problem with the polar- 
ized distributions is restricted to the singlet part of the structure function. 
The same considerations can be drawn from Fig. |, where the polarized d 
valence distributions in LO and NLO are shown together with the data at 
g2 = 10 GeV\ 

Finally, in Fig. ^ we show the difference between the polarized structure 
functions of the proton and neutron (for which the integral is qa)- Only 
the theoretical values are presented. According to (|^), in LO the moments of 
this quantity are just the moments of the parton distribution Au{x) — Ad{x) 
but in NLO we have to include, in addition to the NLO evolution of the 
parton distribution, the second order effects in the form of the corresponding 
Wilson coefficient (listed in the Appendix). The Wilson coefficient leads 
to a correction in the integral over the NLO curve in Fig. ^ of the form 
1 — as{Q'^)/7T. In practice, this reduces the integral from 1.32 in LO to 1.24 
in NLO (as(Q^) was calculated using the asymptotic approximation for equ. 
(D and it has the value 0.191 at = iQGeV^). 



5 Conclusions 

It is clear from the close agreement of the various quark distributions calcu- 
lated in LO and NLO that the excellent results obtained in earlier studies 
based on the bag model were reliable. That the di-quark masses required 
(700 and 900 MeV) were so close to the values expected in the model (600 
and 800 MeV - c.f. sect. 4) is very reassuring. Furthermore, the strong 
coupling constant in NLO is small enough that one can be quite confident in 
the convergence of the calculation. 
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It is important to be clear what has and has not been achieved. As most 
quark models (including the bag) are not derived from QCD one cannot 
calculate structure functions unambiguously at NLO. There is an inevitable 
limit to the accuracy of any quark model. What the present work has shown 
is that once the parameters of the model (e.g. /i^, R, Mg, My) are fixed by 
the unpolarized valence distribution in either LO and NLO, the predictions 
for the spin dependent valence quark distributions are determined relatively 
unambiguously. This is an important result as it does give one some con- 
fidence in our ability to investigate problems like those mentioned in the 
introduction (nuclear EMC effect, spin crisis, etc.) within quark models. 

As an immediate application we plan to estimate the effect of the pion 
cloud of the nucleon on the spin and flavour dependence of the bag structure 
functions p6 . 



We thank to A. Schreiber for providing the program for leading order evo- 
lution which we extended to next-to-leading order. F.M.S wishes to thank 
to W. Melnitchouk for helpful discussions. This work was supported by the 
Australian Research Council and by CAPES (Brazil). 

A Appendix: Beta function, NS anomalous 
dimension and Wilson coeficient 

We here list the coefficients appearing in the expansions (0)-(^. They can 



be found in refs. |28, 29, 31 



P. - H^^^, (A.1) 

A = yC^l - 4 (^^Ca + Cf) TnNj, (A.2) 

^(o)n ^ f-3 - — ^ + 4S,in)) , (A.3) 
V n{n + l) I 
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(l)n 



+64^(n) + 24S2{n) - 



3-8.^(1 
„ 2n+ 1 



3n3 + n2-l 2n2 + 2n + ll 



(n + i: 



+ 52 (n) 


52 






f 160 


2 





n{n + 1) 
32 
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- 165i(r2)52(^) 

151n^ + 263^3 + 97^2 + 3n + 9 ■ 
9ra3(n + 1)3 
lln2 + 5n+ -3l 



/ ISOq^ n , 32,. n , 4 , ^^ lln^ + 5n+- 
j- — + y^^(") + 3 + 9n^(n + l)^ 



1 

-9 + - 



— + 35i(n) -452(n) - 



nin + 1 



(A.5) 



+2{Sl{n) + S^in)) 

For the case of 5'f/(3)c, C^? = 
1 or 2 for analytic continuation of odd or even momenta respectively. The 
functions S{n) are given by: 



Ca = 3 and 



The parameter r] is 



SAn) 



n ^ 



.l)v 



j=i J 



n 



(A.6) 

(A.7) 
(A.8) 



We use the inverse Laplace transformation of the IMSL library to make 
the inversion of expression (|^) for the moments and extract the parton distri- 
butions. Once the moments are defined only for n even or odd, it is necessary 
to extend the validity of the equations for the whole interval of n. This is 
done through analytic continuation in which case the values of the anomalous 
dimensions and Wilson coefficients starting from even or odd n coincide. For 
this purpose, it is better to work with the forms of the S{n) functions given 



in Appendix A of the 1990 paper of Gliick et al. [22 . 



(A.4) 
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Figure 1: Total valence distribution in the bag compared with the MRS 
23] parametrization of the data in the MS scheme. The quark distributions 
are evolved in leading order QCD. 

Figure 2: Total valence distribution in the bag compared with the MRS 
parametrization of the data in the MS scheme. The quark distributions 
are evolved in next-to-leading order QCD. 

Figure 3: Preliminary polarized valence data for the up quark distribution 
in the proton compared with bag model predictions in leading order and 
next-to-leading order. 

Figure 4: Preliminary polarized valence data for the down quark distri- 
bution in the proton p6| compared with bag model predictions in leading 
order and next-to-leading order. 

Figure 5: Bag model prediction for the x dependence of the Bjorken sum 
rule in leading order and next-to-leading order QCD evolution. 
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Figure 1: 
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Figure 3: 
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